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We study the detection of extragalactic point sources in two-dimensional flat simu- 
lations for all the frequencies of the forthcoming ESA's Planck mission. In this work we 
have used the most recent available templates of the microwave sky: as for the diffuse 
' Galactic components and the Sunyaev-Zel'dovich clusters we have used the "Plank 

, Reference Sky Model" ; as for the extragalactic point sources, our simulations - which 

comprise all the source populations relevant in this frequency interval - are based on 
up-to-date cosmological evolution models for sources. To consistently compare the ca- 
pabilities of different filters for the compilation of the - hopefully - most complete blind 
Q ' catalogue of point sources, we have obtained three catalogues by filtering the simu- 

^ I lated sky maps with: the Matched Filter (MF), the Mexican Hat Wavelet (MHWl) 

C/3 . and the Mexican Hat Wavelet 2 (MHW2), the first two members of the MHW Family. 

' For the nine Planck frequencies we show the number of real and spurious detections 

and the percentage of spurious detections at different flux detection limits as well as 
the completeness level of the catalogues and the average errors in the estimation of 
the flux density of detected sources. Allowing a 5% os spurious detections, we obtain 
the following number of detections by filtering with the MHW2 an area equivalent to 
half of the sky: 580 (30 GHz), 342 (44 GHz), 341 (70 GHz), 730 (100 GHz), 1130 (143 
GHz), 1233 (217 GHz), 990 (353 GHz), 1025 (545 GHz) and 3183 (857 GHz). Our 
current results indicate that the MF and the MHW2 yield similar results, whereas 
the MHWl performs worse in some cases and especially at very low fiuxes. This is a 
relevant result, because we are able to obtain comparable results with the well known 
Matched Filter and with this specific wavelet, the MHW2, which is much easier to 
implement and use. 

Key viTords: filters: wavelets; point sources: catalogues, detections 
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1 INTRODUCTION 

In tw o years from now, the ESA's Planck^ satellite llTauberl 
|20o3) will inaugurate a new era in the studies of the Cos- 
mic Microwave Background (CMB) radiation. Planck will 
observe the microwave sky with unprecedented angular res- 
olution and sensitivity in nine frequency channels ranging 
from 30 to 857 GHz. Besides CMB anisotropies, Planck will 
also yield all-sky maps of all the major sources of microwave 

" E-mail: caniego@ifca.unican.es 

^ http:/ /www.rssd.esa.int/index.php?project=Planck 



emission, including a large number of extragalactic sources 
that have not yet been observed at these frequencies. The 
study of these sources of microwave emission -often referred 
to as "contaminants" or "foregrounds"- is twofold: on one 
hand, it is necessary to remove them in order to have the 
cleanest possible view of the CMB and, on the other hand, a 
better knowledge of the different foregrounds is a scientific 
goal in itself. As the launch date approaches, a big efi^ort 
is being made to develop and to test state-of-the-art data 
processing techniques that will optimise the scientific ex- 
ploitation of the forthcoming Planck data. 

The case of extragalactic foregrounds deserves to be 
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considered in detail. Radio and infrared selected galaxies 
will be seen by Planck as point-like objects due to their very 
small projected angular size as compared to the experiment 
resolution (see Table 1). Hence, they are usually referred to 
as extragalactic point sources or just as point sources (as 
we do hereafter). Point sources are expected to be a ma- 
jor contaminant for P l anck at multipoles I ^ 500 — 1000 
(Toffolatti ct al. 1998; Dc Zotti et"al] Il999l: iHobson et all 
Ll999; Tceinar k et al. 20QQ) and it is necessary to detect and 
remove as many of them as possible in order to clean CMB 
maps. As a very important astrophysics by-product, the de- 
tection process will yield all-sky catalogues of extragalactic 
point sources in a frequency interval in which they are lack- 
ing and that will prove very useful to constrain models of 
galaxy formation and evolution. 

Unfortunately, many of the component separation tech- 
niques that are generally used to separate diffuse Galactic 
foregrounds are not well suited to deal with point sources. 
This is mainly due to the fact that each galaxy is an inde- 
pendent source, different, in principle, to any other source 
in the sky. Albeit average energy spectra can be defined 
for different source populations, the spectral emission law is 
different for each galaxy and it makes impossible to apply 
separation methods which exploits a single energy spectrum 
for this foreground component. 

If the lack of knowledge on their spectral emission prop- 
erties can be troublesome, their projected angular shape 
is basically the same for all of them: a Dirac's 5-like re- 
sponse convolved with the instrument beam response func- 
tion. Thus, techniques that take into account the specific 
angular shape of point sources, such as wavelets and band- 
pass filters, are particularly useful for detecting them. In 
the last few years a number of techniques have been pro- 
posed for the specific case of point source detection in CMB 
maps , including th e Mexican Ifat Wavelet (Cavon ct al. 
I2OOOI: IVielva et aP Eo Ola b, 2003), the M atched fiher 
(MF) iTeemark fc de Oliveira-Costal Il99^. the Scale- 
Adaptive filter "^^^jll ^erean^^^Martfnez-GonzalejlioOlt 
[Herranz et al." '2002a"b'), the Adaptive Top Hat filter 
( Chi ang et al. 2002), the Bipar ametric Scale- Adaptive fil- 
ter ilLopez-Caniego et alJ l2005^ an d the recently intro- 
duced Mexican Hat Wavelet Family iConzalez-Nuevo et alJ 
I2OO6I) . All the previously mentioned methods belong to 
the class of wavelet and linear filter techniques. Addi- 
tionally, non-linear techniqu es such as Bayesian detection 
iHobson fc McLachlanI |20o3) have been successfully ap- 
plied to point source detection in the CMB context, but 
since the use of those techniques implies a totally differ- 
ent methodological approach we will focus on the above 
mentioned filters and wavelets. Whereas some works have 
made an effort to compare the performances of some of 
them both theoretically and by usin g almost-ideal si mu- 
lations iBarreirp ct al. 20QiJ Vio. An drcani. & Wamst e'keJ 
I2OO4I: lL6pez-Caniego et alJl2005ir an attempt to compare 
the existent techniques under (almost) "real life" conditions 
for the future Planck mission has not been done, yet. 

In this work we intend to reproduce the conditions of 
a blind point source survey as it will be carried out by 
Planck. We will use the Planck Reference Sky and the nomi- 
nal Planck instrument characteristics and goal performance 
to simulate realistic sky emission as it should be observed 
at the nine frequencies covered by the satellite. Using such 



realistic simulations we will compare the performance of dif- 
ferent filters. The ultimate goal is to decide which is the best 
tool of choice for the Planck case. The 'goodness' of a filter 
will be evaluated according to the following criteria: 

(i) The filter must be well suited to conduct a blind sur- 
vey, that is, it must be able to work well with a minimum 
number of a priori assumptions about the data. 

(ii) Besides, it must be robust against the effect of possi- 
ble systematics. 

(iii) It must yield, after the detection process, a high num- 
ber of positive detections. 

(iv) It must yield, after the detection process, a low num- 
ber of false detections (not higher that, let us say, 5% of 
the integral, or total, number of detections above the corre- 
sponding detection threshold). 

(v) Moreover, additional factors, such as the flux detec- 
tion limit of the output catalogue, its completeness and the 
accuracy with which the positions and the flux densities of 
the sources are estimated, will be taken into account. 

We would like to remind that previous criteria are sim- 
ilar to those required for the Early Release Compact Source 
Catalogue, ERCSC, a blind catalogue of point sources that 
is one of the objectives of the Planck collaboration. In the 
compilation of this kind of catalogue, factors such as quick- 
ness and accuracy of the estimation of the fiux of the sources 
have the priority over a high absolute number of detections. 
Keeping this in mind, criterion | (i) | eliminates from the com- 
petition tools such as the Biparametric Scale- Adaptive filter 
that, even if they are potentially very powerful, require a de- 
tailed knowledge of the probability distribution of the fore- 
grounds^ . The Adaptive Top Hat filter is known to produce 
strong ringing artifacts around the sources which can lead 
to a high number of false detections -against criterion |(iv)| -, 
in particular in the vicinity of bright sources. The Scale- 
Adaptive filte r seems to perform si milarly or slightly worse 
than the MF jHerranz et al.ll2002ij) . Therefore, in this work 
we will focus on the comparison of two tools: the Mexican 
Hat Wavelet Family (MHWF) , which includes the standard 
Mexican Hat Wavelet, and the Matched Filter^. 

At first glance, the MF should be the obvious winner in 
the comparison. By definition, it is the best linear operator 
that can be applied to the data in order to maximise the 
signal to noise ratio of a signal with a known profile embed- 
ded in additive noise. But, in practise, the use of MFs does 
not lack of subtleties that must be considered here. 

In Fourier space the MF is proportional to the inverse of 
the power spectrum of the noise. That means that the power 
spectrum of the noise must either be known a priori or be 
estimated from the data in order to construct the filter. If 
point sources are scarce, the power spectrum of the noise can 

^ Such tools could be very useful for exhaustive data mining the 
sky down to ver y low flux limits. 

^ In particular, iGonzalez-Nuevo et a l. ( 2006) have shown that 
the second member of the Mexican Hat Wavelet Family (to be 
introduced in section |2J is the one that gives the best results for 
the case of the Planck Low Frequency Instrument (LFI) chan- 
nels, and therefore we will limit the discussion to this wavelet, 
comparing it with the standard Mexican Hat (the first member 
of the family), that has been widely used in the literature with 
good results, and the matched filter. 
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be reasonably approximated by the power spectrum of the 
observed data, that is easy to estimate by means of any of 
the power spectra toolboxes available in scientific softwares. 
But any estimated power spectrum, as good as it may be, is 
just a good guess of the real thing. This leads to a number 
of problems: 

• Firstly, it is necessary to estimate the value of the power 
spectrum for all the Fourier modes present in the image. 
This implies the estimation of several hundreds of numbers 
for a typical CMB image. For the typical image size of the 
sky patches we work with, each Fourier mode must be esti- 
mated from a small number of data samples. Therefore, the 
estimated power spectrum is noisy, especially at low Fourier 
modes. On the contrary, to use the various members of the 
Mexican Hat Wavelet Family it is only necessary to deter- 
mine one single number, the optimal scale of the wavelet 
and, thus, it is much less sensitive to noise estimation. 

• Since the estimated power spectrum is noisy, if it is di- 
rectly used to construct the MF it very often happens that 
the resulting filter is so full of discontinuities and 'jumps' in 
Fourier space that strong ringing effects appear in the fil- 
tered image. Therefore, it is necessary to smooth the power 
spectrum before constructing the filter. The different possi- 
ble choices used in the smoothing procedure introduce some 
degree of arbitrarity in the use of MFs: one can, for exam- 
ple, apply some binning and interpolation scheme, or use 
polynomial fitting to the power spectrum, etc. 

• An additional problem appears when it is not possible 
to properly estimate some Fourier modes. This is the case 
when the image which has to be filtered is not complete (for 
example, if a mask is applied to the data in order to cut 
bad pixels, or in areas of the sky where there is incomplete 
coverage by the instrument). In that case the missing modes 
must be somehow guessed in order to build the MF. This 
problem is much less relevant for the case of wavelets. 

• Moreover, if we make considerations in the sphere we 
will have to deal with important problems when using a 
Matched filter as compared with a wavelet. These problems 
arise from the fact that the foregrounds are very different in 
different regions of the sky and therefore we need to use the 
ap propriate filt e r for every region. The approach followed 
bv lVielva etai] tOO^ using the Spherical MHW (SMHW) 
was to divide the sky in a number N of regions and ob- 
tain the optimal scale for all of them. Then they determined 
that many of these scales gave similar results and that they 
could be divided in a small number of different groups. This 
allowed them to construct just a few filters in the sphere 
with their corresponding scales and filter the maps a small 
number of times instead of N. This is important because de- 
pending on the resolution, the filtering process may need a 
lot of CPU time. The problem with the MF is that instead 
of calculating N optimal scales we would need to calculate N 
filters, calculating the power spectra from the N regions, and 
filtering the maps N times, once for every filter. This pro- 
cess will require enormous amounts of CPU time, especially 
when dealing with high resolution images, and therefore will 
make it unfeasible in practise. 

Therefore, and as any experienced practitioner perfectly 
knows, the use of the MF is not the same as filtering by the 
inverse of the squared Fourier transform of the data. On 
the contrary, it requires a non negligible effort of handmade 



tuning that is not free from arbitrarities. Besides, all the 
previous effects lead to an unavoidable degradation of the 
performance of the MF under realistic conditions. 

Furthermore, all the theoretical superiority of the MF 
with respect to other linear filters comes from the fact that 
it maximises the signal to noise ratio, that is, it minimises 
the variance of the filtered noise. But if the noise is not 
Gaussian, as it is the case in CMB maps due to the emis- 
sion contributed by Galactic foregrounds, the fact that the 
variance is minimum does not guarantee that the number of 
false detections be minimum. There may be outliers that are 
not removed. In that case, it is not by any means clear that 
the MF should be better than any other filter. 

Taking all the previous points into account, the com- 
parison between the MF and the wavelets is still necessary. 
The performance of wavelets will degrade as well when go- 
ing from ideal to realistic conditions. Wavelets, however, are 
much less sensitive than MFs to the problems described in 
the paragraphs above. Therefore, it is expected that the per- 
formance degradation will be less severe. If we can find a 
wavelet that performs nearly as well as the MF, but with- 
out having to resort to handmade tuning, we will have a 
detection tool that is as good as the MF regarding crite- 
ria |(iii)| |(iv)| and |(v)| but is better regarding criterion |(i)| 
and |(ii)| Such a wavelet would be preferable to MF for the 
compilation of a Planck blind point source catalogue. 

In Section|^we briefiy review the tools to be used in this 
paper: the MF and the wavelets belonging to the Mexican 
Hat Wavelet Family. In Section |3 we describe the realistic 
Planck simulations we use. The results are summarised in 
Section 01 In Section |S] we discuss some additional issues 
regarding point source detection in microwave satellite mis- 
sions. Finally, we describe our main conclusions in Section|S| 



2 METHODOLOGY 

In this Section we consider the detection of compact (point- 
like) sources where A is the intrinsic fiux density of the 
source and s{x) = At{x) is the observed fiux, filtered by 
a circular Gaussian profile t{x) = exp(— a;^/2f?^), x = \x\, 
whose Fourier transform is r{q) = exp{ — {qR)^ /2), q = 
\q\. Such profile represents a commonly used useful approx- 
imation to the real profile of the beam response function, 
which can be more complicated and not analytically simple 
to describe. Of course, the MF can be easily designed for a 
non-axysimmetric beam, but also elliptical extension of the 
MHWF can be done (for in s tance , the elliptical MHW h as 
been used in lMcEwen et al.l (l2005l) and lCruz et"al] J200el) '). 
In general, the elliptical approximation is good enough for 
the typical beam shapes of the CMB experiments. In any 
case, dealing with non-axysimmetric tools implies increas- 
ing the CPU time, since different orie ntation must be con - 
sidered. Fortunately, as it was shown bv lVielva et alJ i2003l) . 
data with non-axysimmetric but well known main beams can 
be still analysed with isotropic filters by characterising the 
bias introduced for the particular beam shape. 

Given the above, the filtered sources in the map ap- 
pear as circularly symmetric objects and, with the assump- 
tion that the background is statistically homogeneous and 
isotropic, the most natural thing to do - for detecting point- 
like sources - is the application of circularly symmetric fil- 
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ters to the map. More specifically, let us consider a 2D filter 
'^{x; R, b), where R and b define a scaling and a translation 
respectively, then 



(1) 



If we filter our field f{x) with '^'{x; R,b), the filtered map 
will be 



(i?,6) = / dx f(x)^{x-R,b). 



(2) 



The filter is normalised such that the flux of the source in 
the position of the source {b — 0) is the same after flltering: 



dxT{x)^{x;R,0) = 1. 



(3) 



For the filtered map the moment of order-n is defined as 



2tv 



(4) 



where P{q) is the power spectrum of the unfiltered map. 
Note that the zeroth-order moment, ctq, is the dispersion of 
the filtered map. 

As a result, the source profile described above is char- 
acterised by the "natural scale" R. This parameter appears 
in all the filters that we will consider in the following. For 
a given "optimal scale", _Ro, the amplification (A/ao) of the 
source fiux has a maximum value. By working at this partic- 
ular scale - at which a filter maximise the fiux of a compact 
source with respect to the average surrounding background 
value - it is possible to significantly improve the perfor- 
manc e of the chosen filter in terms of detecting compact ob- 
jects JVielva et al.ll2001al l2003t iLopez-Canieeo et ahlliool 



2.1 The Matched Filter (MF) 

The MF can be obtained by introducing a circularly- 
symmetric filter, '^{x; R,b), and imposing the following two 
conditions: (1) (w(7io,0)) — s(0) = A, i. e. w{R,0) is an 
unbiased estimator of the flux density of the source; (2) the 
variance of w{R,b) has a minimum at the scale -Ro, i- e. it 
is an efficient estimator 



tpMF 
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As mentioned in the previous section, the practical imple- 
mentation of the MF requires some additional work. In equa- 
tion ^ the source profile T{q) is known, but the power spec- 
trum P{q) must be estimated from the data. Moreover, the 
normalisation a must be calculated by integrating P{q) . The 
estimation of the power spectrum from the data can be easily 
done, but the estimates are necessarily noisy. In our case we 
use flat patches with sizes varying from 128 x 128 to 512 x 512 
pixels. For small patches like the ones we use, we obtain a 
extremely noisy power spectrum and we need to smooth it 
before it can be used. We have found that a smoothing pro- 
cedure that consists in binning the power spectrum modes 
and performing linear interpolation for the power spectrum 
values lying between the bin centres works well. Our tests 
suggest that the results vary slightly for any number of bins 



between 25 and 50, but the performance degrades quickly 
outside this interval. As for the results discussed in this pa- 
per, we bin the Fourier modes in 40 uniformly spaced bins. 
Even so, the filter that we construct is noisy, as seen in Fi- 
gure In this figure we have compared the three chosen 
filters at each Planck frequency and for one of the simula- 
tions: the MF and the first two members of the MHWF (that 
will be introduced in the next section). It is easy to appre- 
ciate that the resulting MF is very noisy if compared to the 
smooth curves of the two wavelets, that only depend on the 
optimal scale, a quantity which can be easily estimated (see 
below). 



2.2 The Mexican Hat Wavelet Family 

The Mexi can Hat Wavelet Famil y (MH WF) has been dis- 
cussed by iGonzalez-Nuevo et alj (l2006f) and it is the ex- 
tension of the MHW, obtained by applying iteratively the 
Laplacian operator to the Gaussian function. 
The Fourier transform of ipn{x) is 



Tpniq) 



dx xJo{qx)il}n{x), 



(6) 



where q is the wave number, g = |q| and Jo is the Bessel 
function of the first kind. Any member of this family can be 
written in Fourier space as 

^„(Q)cxg^"e-^. (7) 
The expression in real space for these wavelets is 

iP„(x) oc A-^ifix), (8) 

where tp is the 2D Gaussian ip{x) — -^e'^ ^"^ ■ 

Note that il}\(x) is the usual MHW (hereinafter 
MHWl), a wavelet that has been extensively tested in the 
field of point detection with excellent results. Note that we 
caU MHWn the member of t he MHWF with index n . 

The goal in the work bv iGonzalez-Nuevo et af] (l200d) 
was to compare the performance of some members of the 
MHWF and the main conclusion was that the second level 
wavelet ^i(x) did perform better than ^\(x) for the channels 
of the low frequency instrument of Planck. 

In this work we are going to concentrate on these two 
wavelets, i/'i(a::) and xp2{x), and compare them with the 
Matched filter, as described in the previous subsection, for 
both the LFI and HFI channels. 

Now let us consider a field f{x) on the plane R^, where 
X is an arbitrary point. One can define the wavelet coefficient 
at scale R at the point b in the form given by equations Q 
and @. 

The wavelet coefficients Wn{b, R) for the members of the 
MHWF can be obtained in the following form 



Wnib, R) 



dqe "''^ f(q)\pn{qR), 



(9) 



assuming the appropriate differential and boundary condi- 
tions for the field this expression can be rewritten as 

w4b,R)<x J dx[A"f{x)]^(^^^^y (10) 
and the wavelet coefficient at point b can be interpreted as 
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Figure 1. In each panel we compare the three chosen filters in Fourier space as they result at each Planck frequency for one simulation. 
The solid line corresponds to the MHWl, the dashed line to the MHW2 and the dot-dashed one to the MF. The MHWl and MHW2 
operate at their optimal scale, a scale that yields the maximum average amplification of the source flux. The plots show that the MF is 
noisy as compared to the two adopted wavelets. As a remarkable result, we want to point out that the MF and the MHW2 have very 
similar behaviours in most of the plots. 



the filtering by a Gaussian window of tlie invariant (2n)*''- 
order differences of tfie field /. 



3 SIMULATIONS 

We want to compare the capabilities of the MHWl, the 
MHW2 and of the MF when dealing with the detection of 
compact sources in realistic CMB maps. 

To test these tools in the most realistic conditions so 
far, we used the latest "Planck Reference Sky Model" pro- 
vided by the Working Group 2 (WG2, "Component sep- 
aration") of the Planck Consortium. All the available as- 
trophysics foregrounds compon ents are provided i n all-sky 
maps m the HEALPix format llGorski et alJl2005ll and in- 
clude the following components: 

• Thermal dust emission: has been created using a combi- 



nation of two grey bodies with mean emissivity parameters 
Qi — 1.67 and 02 = 2.70 and mean temperatures Ti = 9AK 
and T2 = 16. 2A" based on the lFinkbeiner. Davis. &: Schlegel 
41999) model. 

• Synchrotron emission : is a cleaned ver sion of the 408 
GHz iHaslam et all lll982h map made by iGiardino et alJ 

• Free-Free radiation: it adopts the free-free model of 
ijDickinson, Dayies, & Davis (_2Q03) which ha s been obtaine d 
by differe nt Ha surveys (e. g. WHAM jHaffned Il999f) . 
SHASSA JCaustad et al.ll200lD ). 

• S-Z Clusters : fol l ow-up of the works by 
IColafr anccsco et al.l lll99'?r) : lDe Zotti et al.l (Il99l|200l) and 
by P. Mazzotta. 

As for point sources (not contained in the "Planck Re- 
ference Sky model", yet) we have used the simulation soft- 
ware (EPSS-2D) presented bv lConzalez-Nuevo etail (l2005l) 
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Figure 2. In this figure we sliow one of tlie simulations at 70 
GHz. In the upper panel we show a simulation of point sources 
after they have been convolved with the beam of the instrument. 
In the middle and lower panels we have the the combination of the 
different components before and after filtering with the MHW2. 



which is a very fast and efficient tool for simulating maps of 
extragalactic point sources. For this first comparison among 
the capabilities of different filters for source detection, we 
have adopted a simple Poisson distribution of point sources 
in the sky, not considering the clustering properties of the 
different source populations. 

The input number counts are the ones foreseen by 
the lDe Zotti et all l |2005|) evolution model for ra dio-selected 
sources and by the lOranato et al.l (|200ll l2004 ) model for 
high-z proto-spheroidal galaxies. These two models have 
proven very successful in fitting the observed source num- 
ber counts recently obtained by the VSA, ATCA, WMAP 
and CBI surveys (at ~15-30 GHz), and by the SCUBA and 
MAMBO surveys (in the far-IR/sub-mm domain), respec- 
tively. As for late-type far-IR gal axies we followed the same 
phenomenological approach as in iNegreHo et al.l J2004) . 

At each frequency, 145 patches of 12.8x12.8 square de- 
grees have been simulated, covering a total effective area 
of 27r sr (half of the sky). Each patch has contribution 
from: a CMB realisation following the Concordance Model 
iSpergel et al.i r2003'l. a point source realisation given by the 
model described above and the sum of all the Galactic com- 
ponents. Galactic foregrounds have been selected from 12 
different regions of the sky (at Galactic latitude b > 30°). 
These 12 regions are representatives (in terms of its disper- 
sion) of the Galactic emission outside the above mentioned 
Galactic cut. Each one of these 12 Galactic simulations has 
been used accordantly to their represe ntativeness. This is a 
similar procedure to the one followed in lVielva et"ai] (l2001a^ 
(in that case 10 Galactic patches were considered). 

Once the CMB, the point sources and the Galactic com- 
ponents have been added, the resultant map is convolved 
using the instrumental Planck specifications (see Table 0. 
Finally a Gaussian white noise realisation is added, with a 
dispersion per pixel according to the Planck specifications 
(see Table 1). 

Each image is filtered with the MF, the MHWl and the 
MHW2. To avoid spurious artifacts due to border effects, 
five pixels on the edge of the image have not been consid- 
ered. The MHWl and MHW2 operates at their respective 
optimal scales. This optimal scale can be easily obtained by 
filtering the image with a series of scales and by finding the 
particular scale that produces an image with maximum am- 
plification. Note that this procedure does not affect the fiux 
of the sources, since the filters are normalised in such a way, 
eq. (3), that the flux of the source is always preserved. 

We then compare the list of source candidates (our re- 
covered catalogue) with the reference (input) source cat- 
alogue, by trying to match the detected sources with the 
simulated ones. A source candidate is considered as a pos- 
itive detection if it is found within a radius r ^ FWHM/2 
from the corresponding source in the reference catalogue. 
Only those objects in the candidate catalogue with a pos- 
itive matching and with an error in the estimate of their 
fiux density ^ 100% (with respect to the fiux of the simu- 
lated source) are considered as real (positive) detections. Ev- 
ery other candidate source with flux above a fixed detection 
threshold in the filtered map (see Figures ) is considered as 
a false or spurious detection and it has not been included in 
the final catalogue, corresponding to that particular Planck 
channel. 
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Freq. Im. Size FWHM Fix. Size Unoise 
(GHz) (pixels) (arcmin) (arcmin) (10~^) 



857 


512 


5.0 


1.5 


6700.0 


545 


512 


5.0 


1.5 


147.0 


353 


512 


5.0 


1.5 


14.7 


217 


512 


5.0 


1.5 


4.8 


143 


512 


7.1 


1.5 


2.2 


100 


256 


9.5 


3.0 


2.5 


70 


256 


14.0 


3.0 


4.7 


44 


128 


24.0 


6.0 


2.7 


30 


128 


33.0 


6.0 


2.0 



Table 1 . Instrument performance goals for all the Planck satellite 
channels. The antenna is assumed to be a circular Gaussian one 
and the instrumental noise level is in units of AT/T per FWHM 
(e.g in a square whose side is the FWHM of the beam). 



4 RESULTS AND DISCUSSION 

In this Section we compare the performances of the filters 
described and discussed in the previous sections (MHWl, 
MHW2, MF). 

In Figuresrnto llll we present a number of relevant plots 
for all the Planck frequency channels. In the upper panels of 
each Figure we show the integral (or total) number of real 
and of spurious detections (left and right panels, respec- 
tively) above different flux detection limits in half of the sky 
(see captions for more details). It is easy to appreciate that 
the three filters detect approximately the same number of 
objects in all the Planck channels. As regards the spurious, 
or false, detections, the MF and the MHW2 yield very simi- 
lar results whereas the MHWl detects a significantly higher 
number of spurious sources in comparison with the other 
two filters, in particular at frequencies ^ 100 GHz and for 
faint flux detection limits. 

In the lower left panels of Figures 01 to 1 1 II we show the 
percentage of spurious detections with respect to the integral 
number of detections. These plots are very relevant ones 
because they compare the relative performance of the filters 
when taking into account both the number of real and of 
spurious detections. By these plots we can also establish the 
flux detection limit above which it is possible to reach a given 
percentage of spurious, e.g. 5%. It can also be noted that the 
behavior of the MF and the MHW2 is always comparable, 
whereas the MHWl give worse results than these two filters 
at higher frequencies. At the lowest frequencies of the Planck 
experiment the MHWl detects a relative higher number of 
spurious sources, but only when going down to faint flux 
detection thresholds. 

In the lower right panels we also plot the completeness 
level of the catalogue above different flux detection thresh- 
olds. Each flux threshold represents the flux at which we 
have detected a given percentage (e.g., 95%, 100%) of all 
the simulated sources in the map. At low and intermediate 
Planck frequencies we find that our catalogues are 95% com- 
plete (the 95% of the simulated sources are detected above 
this flux limit) above ~ (300 — 350) mjy whereas at the high- 
est Planck frequency, 857 GHz, the 95% completeness level 
is obtained only for fluxes ^ 550 mJy. In Table |21 we have 
summarised these results for each Planck frequency chan- 
nel and we have also included the 100% completeness level. 



the flux limit at which all the simulated point sources are 
detected. 

In Figure [T^ we have plotted, for each Planck channel, 
the number of spurious detections against the number of real 
detections. This kind of representation, together with the 
plot of the percentage of spurious detections (in the lower 
left panels of Figures |21 to IIH , is very important because 
it shows the global performance of the fllters. It is clear 
that, for a given number of spurious detections, the filter 
which is able to obtain the corresponding highest number 
of real detected sources should be considered the best one. 
In this scheme, a filter performs better than another if its 
curve tends to lie above the other ones. In the upper panels 
of the Figure (i.e., at low Planck frequencies) we observe 
that the MF and the MHW2 behave similarly and better 
than MHWl, whose curve is always below the other two. 
At intermediate frequencies the three filters here discussed 
perform similarly, and at 217 GHz and 353 GHz the MHWl 
curve is very close to those representing the MF and the 
MHW2. At 857 GHz, the curve representative of the MHWl 
is again below the other two ones, but only at fluxes above 
~ 650 mJy. 

In Figure 1131 we represent the average error in the 
estimate of the flux of the detected sources above dif- 
ferent flux detection limits, and this error is defined as 
< {Aest — Areai) / Are.ai >■ This average error can take up 
negative values and instead of including the standard devi- 
ation we have plotted the true 68% error bars taking into 
account the probability distribution. 

As expected, the error bars decrease towards higher 
fluxes. In Figure 1131 it can be seen that there is a positive 
bias in the estimation of the flux. Such bias is negligible for 
high fluxes and increases towards lower fluxes. The bias is 
due to the well-known selection effect that makes more likely 
to detect a source when it is located on a positive fluctuation 
of the background, leading to a greater number of overesti- 
mated fluxes than the number of subestimated fluxes. This 
can also be appreciated in the non-symmetric shape of the 
observed error bars. 

In Figure riTI we have represented the average distance 
between the position of the detected sources and their ac- 
tual position in the simulated map before filtering. We have 
experienced that in some cases, after filtering, the maximum 
of a detected source has moved to an adjacent pixel. This 
fact would normally not affect our results, when matching 
candidate objects with the reference catalogues, because of 
the way we compare the simulated (input) and recovered 
(detected) sources (allowing a radius of ^ FWHM/2 for 
the proper matching). In any case, in order to quantify this 
effect, we have calculated this distance for every detected 
source in each map. We have found that, even in the worst 
cases and only at very low Planck frequencies (i.e., 30, 44 
GHz), the average distance is always smaller than one pixel. 

As for the MHWl, we have also checked to what extent 
a multi-fit approximation for estimating the source flux (i.e. 
a fit of the wavelet coefficients to four scales 0.5, 1, 2 and 
3 times the optimal scale) can improve the number of real 
dete ctions, as discussed in previo us works on the same sub- 
ject iVielva et alJl2001allbll20o3) . At frequencies above 100 
GHz, the percentage of spurious detections with this wavelet 
decreases to the level of the MHW2, and at 857 GHz they 
are quite similar. On the other hand, at Planck frequencies 
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below 100 GHz (LFI), while we still find a small improve- 
ment in comparison with the results previously plotted, the 
percentage of spurious detections does not reach the level of 
the MHW2, which performs always better. 

As a final comment, it is clear that the multi-fit approx- 
imation only improves the estimation of the source fluxes 
at the highest frequencies and, correspondingly, reduces the 
number of the spurious detections. However, only at 857 
GHz the performance of the MHWl - with this method - 
is similar to that of the MHW2. As for the error in the 
estimated flux of the detected sources, we found that this 
method applied to the MHWl yields estimated errors simi- 
lar to those determined for the other filters. 



5 FURTHER DETAILS ON THE DETECTION 
OF POINT SOURCES 

In this section we briefly discuss some additional issues re- 
lated to point source detection. 

5.1 Confused/ blended extragalactic point sources 

Every time we observe the sky with a pencil- or a synthe- 
sised beam, if two or more point-like sources (i.e., distant 
galaxies or clusters) lay inside the beam area they appear 
as indistinguishable since they "blend" into a single source: 
a single "spot" in the map having the same angular dimen- 
sion (FWHM) of the beam. As a result, only one "source" - 
whose flux is the integrated value of all the individual sources 
in the beam - is ac tually detected. How ever, if source counts 
are not very steep jPe Zotti et alJl99(J) . this integrated flux 
value is mainly determined by the brightest source in the 
beam, whereas the much more numerous faintest sources 
are o nly slightly modifying the total flux an d adding some 
noise iFranceschini et al]ll989l : lBarconj|l992^ . Thus, source 
blending is unlikely to affect the counts of bright sources, 
which are the ones actually detected. 

As an example, at 30 GHz, if we reach the flux detection 
limit at which we have 5% of false detections, ~ 0.4 Jy, we 
have less than 0.007 sources per beam area and, thus, the 
probability that two such sources fall inside the same beam 
is almost negligible. Moreover, at higher Planck frequencies 
the situation is still better. At 100 GHz we find that at the 
5% false detections limit, the average number of sources per 
beam is ~ 0.0014, whereas for the 857 GHz channel that 
number decreases to ~ 0.0008 sources per beam. On the 
other hand, we have to remind that our simulations do not 
take into account the clustering of sources. If source posi- 
tions are correlated in the sky - as it is well known for all the 
source populations observed by Planck - the probability of 
finding two bright galaxies at the same angular separation 
in the sky increases. Therefore, with a low-resolution exper- 
iment, like Planck, we will observe the summed contribution 
of groups (or clumps) of sources within the sa me beam area. 
The p roblem heis been recently discussed bv lNeerello et alJ 
(|2003), who found that, in the presence of strong clustering, 
source confusion can be important even above the canonical 
5(T detection limit, thus modifying source number counts. 
However, this effect is likely to be important only at Planck 
HFI frequencies and for the highly clustered high-redshift 
proto-spheroidal galaxies (see, e.g., the thorough analysis by 



iGonzalez-Nuevo et alJ (120051') '). In any case, this study goes 
beyond the purposes of this paper and it shall be tackled in 
a forthcoming one. 

5.2 Other compact sources 

Other compact but not completely point-like objects will 
be observed by Planck. The most important among these 
objects are galaxy clusters observed thanks to the Sunyaev- 
Zel'dovich (SZ) thermal effect. Detailed simulations of the 
local universe have shown that the SZ signal from low- 
redshift clusters and superclusters should be marginally de - 
tectable by Planck llDolag et alJl2005t [Hansen et alJl2005^ . 
On the other hand, very few SZ clusters show angular scales 
larger than ten arcminutes. Since the scale of these ex- 
tended clusters is significantly different from the scale of the 
point sources (i.e., the beam size), the filtering with either a 
wavelet or a matched filter - both optimised for the scale of 
the point objects - will tend to cancel the SZ signal just in 
the same way as it cancels the rest of the diffuse signal and 
the large scale fluctuations of the background. Therefore, 
large SZ clusters are not expected to affect signiflcantly the 
detection process. 

Most of the SZ clusters, however, shall be observed 
as "unresolved" or point-like (i.e., smaller than the Planck 
beam). In the latter case, they have to be treated in the same 
way as point sources, except for the fact that their fluxes 
are typically much fainter. Even the brightest clusters are 
on the verge of the detectability, given the sensitivities of 
Planck detectors. As already discussed by many authors, it 
shall be necessary to combine the signal measured in various 
frequency channels in order to detect clusters with Planck 
(Hcrranz et al. 2002ci,). When using a channel- by-channel de- 
tection strategy, as the one presented in this work, the im- 
pact of SZ clusters is very small. In fact, ai v < 217 GHz, 
SZ clusters appear as "cold spots" and therefore cannot be 
detected above any positive flux threshold. At 217 GHz, the 
thermal SZ effect is zero and the kinematic SZ effect is well 
below the level of the CMB fluctuations. Above 217 GHz 
there will be some minor effect due to the thermal SZ effect 
that may lead to a small number of false detections. The de- 
tailed study of this kind of contamination shall be addressed 
in a future work. 

As for Galactic compact sources, representing different 
evolutionary stages in stellar evolution - i.e., compact pre- 
stellar cores, young stellar objects and supernova remnants 
- many of them shall be observed by Planck. Current esti- 
mates predict that Planck surveys should detect from many 
hundreds to thousands of these sources, depending on the 
subclass (see Table 5.3 in the Planck BlueBook, "The Sci- 
entific programme", ESA-SCI(2005)1 In particular, one 
of the most relevant classes of Galactic compact sources is 
that of HII regions, which provide relevant information on 
early stages of stellar evolution and on the Galactic spiral 
structure. The thorough analysis made by iPaladini et alJ 
(2003) has demonstrated that the vast majority (> 80%) 
of the 1442 HII regions in their Synthetic Catalogue shall 
be detected by Planck with a very high signal-to-noise ra- 
tio {S/N ^ 100). In any case, all these sources have well 

* http:/ /www.rssd.esa.int/index.php?project=PLANCK 
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Freq(GHz) 


Filter 


Nd5 




S5(mJy) 


SNRs 






S95(mJy) 




Sioo(mJy) 


30 


MHWl 


497 


26 


448 


3.79 


12 


583 


329 


327 


545 




MHW2 


580 


30 


391 


3.75 


12 


595 


324 


432 


447 




MF 


594 


31 


382 


3.63 


10 


594 


324 


432 


447 


44 


MHWl 


284 


16 


505 


4.10 


12 


388 


369 


263 


495 




MHW2 


342 


18 


443 


3.99 


12 


399 


361 


263 


495 




MF 


335 


17 


448 


4.07 


10 


418 


349 


303 


447 


70 


MHWl 


304 


16 


510 


4.59 


8 


445 


339 


276 


495 




MHW2 


341 


18 


473 


4.54 


9 


444 


339 


276 


495 




MF 


344 


18 


467 


4.49 


9 


443 


339 


276 


495 


100 


MHWl 


586 


30 


294 


4.45 


6 


1237 


138 


871 


198 




MHW2 


730 


38 


248 


4.42 


7 


1323 


131 


877 


197 




MF 


690 


36 


252 


4.66 


6 


1325 


131 


877 


197 


143 


MHWl 


914 


48 


240 


4.80 


7 


1713 


116 


1048 


195 




MHW2 


1130 


59 


196 


4.73 


7 


1797 


109 


1049 


195 




MF 


829 


43 


250 


6.25 


2 


1714 


116 


1049 


195 


217 


MHWl 


1153 


61 


195 


4.87 


6 


1756 


121 


1074 


195 




MHW2 


1233 


65 


185 


5.11 


5 


1755 


121 


1074 


195 




MF 


1242 


65 


195 


5.57 


3 


1755 


121 


1074 


195 


353 


MHWl 


887 


47 


314 


5.29 


5 


2253 


147 


1161 


245 




MHW2 


990 


52 


292 


4.86 





2243 


147 


1160 


245 




MF 


986 


52 


292 


4.99 





2247 


147 


1160 


245 


545 


MHWl 


726 


38 


785 


6.33 


8 


2582 


349 


1495 


495 




MHW2 


1025 


54 


647 


5.09 


9 


2568 


349 


1495 


495 




MF 


924 


49 


684 


5.60 


9 


2580 


349 


1495 


495 


857 


MHWl 


1034 


54 


1765 


8.24 


2 


13885 


532 


5194 


745 




MHW2 


3183 


167 


1013 


4.82 


7 


13851 


532 


5201 


745 




MF 


3093 


162 


1020 


4.97 


6 


13914 


532 


5200 


745 



Table 2. In this table we summarise the 5% percentage of spurious detections and the 95%/100% completeness of the catalogues at 
the Planck frequencies. Regarding the percentage of spurious detections, columns three to seven show the number of detected sources, 
Afjjg, the number of spurious sources Ngp^, the corresponding flux detection limit, S5, the signal-to- noise ratio SNR for Ss and the 
average error in the estimation at this flux limit, respectively. We allow for a 5% of spurious detections, where this percentage has been 
calculated as Ngp/{Nj^ + Ngp). Regarding the completeness level, the flux limits at which the catalogues are complete at the 95% and 
100% are shown in columns eigth and ten, respectively, and their corresponding number of detections, N^ig^ and N^-^^^ , in columns nine 
and eleven, respectively. 



known positions in the sky and, moreover, almost all of them 
lay at very low Galactic latitude. Therefore, they shall be 
easily identified and subtracted out. Finally, Planck should 
also observe 397 asteroids and a siRnifican t fraction (50-100 
objec ts) of them with a high S/N ratio iCremonese et all 
|2003)- However, these Solar System bodies are moving ob- 
jects very close to the Earth and, thus, they shall show vari- 
able positions and fluxes, depending from their heliocentric 
and geocentric distances, even during a single Planck scan 
of 360°. 



5.3 Systematics 

As explained in section |3| our simulations can be regarded 
as "state of the art" as for the realism of the simulations of 
the physical components (CMB, Galactic foregrounds and 
extragalactic objects) and we have used the nominal noise 
levels, beams and pixel sizes of the Planck mission. However, 
some systematics of relevance have not been included in the 
simulations. Among them, the ones that can affect our cur- 



rent results are: anisotropic noise, 1// noise, thermal effects 
and non-symmetric beams. 

1/f noise and thermal effects in the electronics of the 
satellite will lead to large scale noise fluctuations across the 
sky. Since the fllters here discussed operate at small angular 
scales, these large scale fluctuations tend to be cancelled 
out by the filters and should not affect the detection of 
point sources. Anisotropic noise will appear due to the non- 
uniform sky coverage of the satellite instruments. Further- 
more, mapmaking algorithms can introduce some degree of 
noise correlation over small angular scales. Non-uniformity 
of the noise will make easier to detect point sources in some 
areas of the sky (those that are better scanned) and more 
difficult to detect in others: on average, the integral number 
of detections should not change very much. Noise correla- 
tion, on the other hand, can be a non negligible effect and it 
must be studied in a future work. In any case, IVielva et alJ 
I2OO3I) have found that their results with the Mexican Hat 
Wavelet did not vary when anisotropic noise pattern and 
1/f noise were considered instead of isotropic white noise. 

We have worked with idealised beams that are circularly 
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symmetric and Gaussian-shaped. In a real experiment this 
will not be the case. The real response of the beam shows 
more complexity, with side lobes, ellipticity and even chang- 
ing the orientation in the sky of the projected shape. As long 
as the beam profile is well known, it is straightforward to 
construct the corresponding MF. Even in the case of non- 
symmetric beams this can be done, by increasing the compu- 
tational complexity of the problem. Therefore, when using 
the MF, real beams can be in prin ciple handled . Rega rding 
the standard Mexican Hat Wavelet. I Vielva et alJ i2003tl have 
tested the influence of realistic asymmetric beams. Those au- 
thors have shown that, although the MHW is an isotropic 
wavelet, it can also be adequate to perform the detection 
of point sources that show a slight Gaussian asymmetry. 
This is precisely the situation for the Planck beams (see the 
Planck BlueBook, ESA-SCI(2005) 1). For a m ore detailed 
discussion, the reader is referred to lVielva et al.. (■2003.'1 . 

5.4 Extension to the sphere 

If we want to extend our work to the sphere, it is important 
to address the problem of very variable foregrounds from 
one region of the sky to the other. Then it is necessary to 
operate locally. From the point of view of implementation, 
there are two approaches to this problem: to work with small 
fiat patches, as we have done in this paper, or to work di- 
rectly on the sphere using the harmonic transform instead o f 
the Fourier transform, as it was done in lVielva et alJ i2003l) . 
The first approach is valid if the patches are small enough 
(no larger than ~ 15 x 15 square degrees), but it implies to 
project the data from the sphere to the plane and to repix- 
elise the samples (from some non-square pixels on the sphere 
such as the ones used by HEALPix to the square fiat pixels 
of the patches). This leads to some deformation, especially 
for the outer parts of the flat patches and for high latitudes 
in the sphere. Therefore, the flat patch approach leads to 
some amount of "systematic" degradation of the data. Be- 
sides, in order to cover all the sky it is necessary to decide 
how to choose the positions of the patches, how much over- 
lapping among the different patches is necessary and how to 
correct the estimated fiux of each source taking into account 
the deformation just mentioned. 

Therefore, in some cases it could be preferable to work 
directly in the sphere. The most simple way to take into 
account the variability of the foregrounds is to use the har- 
monic transform to filter the sphere several times, using each 
time a filter that is adapted to the statistical properties 
of a given area of the sky. This is the approach adopted 
bv IVielva et alJ i2003r) . But the number of different filters 
used should not be very large or the CPU time required 
would be unreasonably long. Fast harmonic transform al- 
gorithms make possible to filter the whole sky in a rela- 
tively short time (Gorski ct al. 2005), but even so the pro- 
cess can take a significant amount of time for high resolution 
maps (for example, performing the harmonic transform and 
calculating the angular power spectrum of a NSIDE=2048 
HEALPix map takes 108 minutes in a 3 GHz Intel Pen- 
tium processor). Thus, it is necessary to reach a compromise 
between acting locally (that may require to filter the data 
many times) and saving CPU time. 

In the case of wavelets it has been successfully shown 
JVielvaet al.ll2003h that it is possible to estimate locally 



the optimal scale of all the regions, group together different 
regions with similar optimal scales and construct a small 
number (around fifteen in the worst case) of filters on the 
sphere for these optimal scales. In their study it was possi- 
ble to reduce the number of times it was necessary to filter 
the maps due to the fact that the functional form of the 
wavelets depend only on one parameter (the optimal scale) 
that varies slowly across large zones of the sky and it is easy 
to group the optimal scales in a small number of bins with- 
out loosing much efficiency. Even so, their algorithm takes 
~ 72 hours for the most CPU time-demanding channel (857 
GHz) in a Compaq HPC320 (Alpha EV68 1-GHz proces- 
sor) and requires 4 GB of RAM memory, a non-negligible 
amount of computer resources. 

On the contrary, the MF functional form depends on the 
power spectrum of the data. The shape of the local power 
spectrum varies wildly from one region of the sky to other, 
and it is not easy to group local power spectra in classes 
to reduce the number of times the sky must be filtered. A 
taxonomic criterion to group power spectrum curves should 
be provided, introducing a new arbitrarily in the process. 
If one wants to skip this arbitrarily, it is necessary to con- 
struct one filter for every region of the sky and to filter the 
sky map a large number of times, especially for high reso- 
lutions. Nowadays, this process requires a huge amount of 
computational resources. This gives another reason to prefer 
wavelets to MF when working in the sphere. 



6 CONCLUSIONS 

We have compared the performance of three filters when 
dealing with the detection of point sources in CMB fiat sky 
maps. These filters are the well known MF, the Mexican 
Hat Wavelet (MHWl) and the recently introduced Mexican 
Hat Wavelet 2 (MHW2). The latter is the second member 
of the so-called Mexican Hat Wavelet Family (the MHWl 
is the first), a group of wavelets obt ained by applying the 
Laplacian operator to the Gaussian iGonzalez-Nuevo et alJ 
2006). This new wavelet is part of an effort to improve 
the MHWl, a tool already exploited by our group that has 
proved very useful in the detectio n of compact s ources and 
of no n-Gaussianity in CMB maps iVielva et al JioOla, 200i 
|2004 . 

We have tested these tools in realistic 2D simulations of 
the microwave sky prepared in the framework of optimising 
the efficiency in separating the various astrophysical com- 
ponents in the forthcoming ESA's Planck Satellite all-sky 
maps. As for the Galactic foregrounds and the S-Z effect 
in clusters of galaxies, we have used the available "Planck 
Reference Sky Model" ; we have adopted the Standard " con- 
cordance" Model for simulating CMB anisotropy maps and, 
as for extragalactic point sources, we have used simulations 
obtained by the software discussed in lGonzalez-Nuevo et alJ 
^2005^ and with the n u mber count models for so urces 
of iDe Zotti et al.l (|2005^ : iGranato et all J200ll. 12004) and 
lNegren^^^Ln2004) . We then applied the three considered 
filters to a sufficient number of fiat patches to cover half of 
the sky (27rsr, 6 > |30°|) (see Section 3). 

We have found that the MHW2 and the MF outperform 
the MHWl in some aspects, especially at the lowest Planck 
frequencies. The three filters yield approximately the same 
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number of real detections, down to the same flux detection 
limit, although the MHWl yields a corresponding higher 
number of spurious detections. Moreover, the MHW2 and 
the MF give comparable results for almost every one of the 
analysed indicators. As shown in Figure it is remarkable 
that, even the estimated shape of the MF tends to that of 
the MHW2, at its optimal scale, for most of the frequencies 
discussed in this work. 

In a previous work JVielva et alJl2003D . a multi-fit ap- 
proximation to estimate the flux density of the sources de- 
tected by the MHWl was shown to be able to improve the 
results (as compared with the direct approach used in this 
work) . We have applied this technique to the simulated maps 
used in this work and we found different results for the LFI 
and the HFI Planck frequency channels. In the first case, 
LFI, the improvement is small and the final results are never 
comparable to the ones obtained with the MHW2. In the sec- 
ond case, this procedure helps to slightly reduce the number 
of spurious sources, except for the highest frequency of HFI 
(857 GHz), where the decrease in the number of spurious is 
significant and the MHWl approaches the results obtained 
with MHW2. 

These results are very important because both wavelets, 
the MHWl and the MHW2, are represented by a known 
analytical function. The only parameter that needs to be 
obtained for each simulation is the "optimal scale", which 
is calculated locally in a very easy way. On the other hand, 
the correct definition of the MF implies a number of steps. 
Firstly, it is necessary to estimate the value of the power 
spectrum for all the Fourier modes present in the image, 
especially for the low modes where the power spectrum is 
noisy. Secondly, the use of such a noisy power spectrum to 
construct the MF often yields a filter with many discontinu- 
ities in Fourier space which, in turn, produces ringing effects 
in the filtered image. Therefore some smoothing in the spec- 
tra needs to be done before constructing the filter, which 
introduces further arbitrariness. Thirdly, sometimes it will 
not be possible to properly estimate some Fourier modes, 
for example when using masks with missing data, and these 
modes will have to be guessed. 

Therefore, the most relevant conclusion of this analysis 
is that the MHW2 can be surely a better choice for the defi- 
nition of a blind source catalogue - like the ERCSC foreseen 
for the future Planck mission - because it gives numbers of 
detected and of spurious sources comparable to the ones 
obtained with the MF but is easier to implement, more ro- 
bust and has much lesser CPU requirements, especially in 
all-sphere applications. 
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Figure 3. 30 GHz. In this figure we compare the results obtained for the three considered filters, the MHWl and MHW2 at their optimal 
scales and the MF. In the upper left figure we show the number of detected sources above a given flux detection limit. In the upper right 
figure we present the number of spurious or false detections above a given flux detection limit. The lower left flgure shows the percentage 
of spurious sources (deflned as the number of spurious divided by the sum of spurious and real detections) above a given flux detection 
limit. The lower right figure shows the completeness level of the catalogues above a given flux detection limit. The solid line corresponds 
to the MHWl, the dashed line corresponds to the MHW2 and the dot-dashed line corresponds to the MF. Similarly, the circle, cross and 
squares correspond to the MHW, MHW2 and the MF respectively. 




Figure 4. The same as in Figure Islbut at 44 GHz. 
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Figure 6. The same as in FigurcElbut at 100 GHz. 
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Figure 8. The same as in Figured but at 217 GHz. 
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Figure 9. The same as in Figured but at 353 GHz. 
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Figure 10. The same as in Figure 1^ but at 545 GHz. 
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Figure 11. The same as in Figure 1^ but at 857 GHz. 
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Figure 12. In each panel we plot the number of spurious detections vs. the number of real detections for all the Planck frequency 
channels. This plot is very useful for comparing the performance of the chosen filters. 
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Figure 13. For each Planck channel we show the average error (68% error bars) in the estimated flux density for all the sources observed 
above different flux detection limits. The error is defined as (< Aest — Areal)/^real >i where Aest is the estimation of the flux in the 
simulated patch after filtering and A^^^i is the flux of the source in the catalogue. For the sake of clarity, we have plotted the errorbars 
with an offset of 10 mjy in the x-axis. 
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Figure 14. In each panel and for different flux detection limits, we sfiow the average distance between the centre coordinates of the 
detected source with respect to the actual coordinates of the corresponding simulated source in the input catalogue. For the sake of 
clarity, we have plotted the errorbars with an offset of 10 mjy in the x-axis. 



